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It is demonstrated that a k-dependence of the hybridization matrix element between /- and 
conduction electrons can give rise to an anisotropic hybridization gap of heavy fermions if the 
filling of electrons corresponds to that of the band insulator. The most interesting case occurs 
when the hybridization vanishes along some symmetry axis of the crystal reflecting a particular 
symmetry of the crystal field. The results of a model calculation are consistent with wide range 
of anomalous properties observed in CeNiSn and its isostructural compounds, the anisotropic 
semiconductor of heavy fermions. In particular, highly sensitive effect of impurity scattering 
on the residual density of states for zero energy excitation and the anisotropic temperature 
dependence of the resistivity are well explained. It is also discussed that a weak semimetallic 
behavior arises through the weak k-dependence of the /-electron self-energy S / (k, 0) . 

KEYWORDS: heavy fermion, Kondo insulator, anisotropic hybridization gap, crystal field, CeNiSn, CeRhSb 



1 



§1. Introduction 

Heavy-fermion systems exhibit a variety of ground states. Among them there is a class, the so- 
called "Kondo insulator" , which has a narrow energy gap of the order (~ WK) at low temperatures. 
This class of compounds exhibits the Kondo effect at high temperatures and an insulating behavior 
at low temperatures .0!) The mechanism of the energy-gap formation at low temperatures has 
been discussed from a variety of view points. These are classified roughly into two categories, i.e., 
k-space and real-space approaches. In the former approach, the origin of insulating behaviors is 
attributed to a hybridyzation gap which is highly renormalized by the strong correlations among 
almost localized /-electrons^ A circumferential evidence supporting this point of view is that all 
the compounds called the "Kondo insulator", except for TmSe, have even number of electrons in 
the unit cell which is a necessary condition for the band insulator. In the latter approach, on the 
other hand, it is attributed to the formation of local bound state of one kind of another, such as 
local singlet due to the strong Kondo effect, the Wigner crystal formation, and so on.B& 

A picture of the renormalized hybridization gap is based on the principle of "adiabatic continu- 
ation"& which was applied first by Landau in the Fermi liquid theory^ and then has turned out 
very useful so far in understanding the low energy properties of wide range of strongly interacting 
systems as demonstrated by Yamada and Yosida for the impurity Anderson model.& While cor- 
rectness of this picture for the "Kondo insulator" has been suggested by the Gutzwiller approach 
for the periodic Anderson model for this decade,!' 0) it was recently shown more vividly on the 
basis of the numerical renormalization group method©, the quantum Monte Carlo calculation!!^ 
and the exact diagonalization method© with the help of the d = oo theory. On the other hand, it 
is reported that a charge gap is different from a spin gap in the d = 1 theory on the basis of the 
exact diagonalization method© This result is against the picture mentioned above. However, we 
believe that this may come from the particularity of 1-dimensionality. 

Among the "Kondo insulators", CeNiSn and its isostructural compounds have attracted much 
attention for almost a decade, because they have exhibited, at least for samples of early stage, 
behaviors of anisotropic semiconductor with highly renormalized energy gap.00'© An anisotropic 
semiconductor means that the energy gap vanishes in a direction in the k-space and the density 
of states (DOS) shows a V-shaped like structure at low energy region. This was first inferred 
from the measurements of the longitudinal NMR relaxation rate at low temperatur MEM) and 
reinforced later by those of the thermodynamic propertiesciUfJ (the specific heat, the magnetic 
susceptibility etc.). In addition, the inelastic neutron scatterings show the existence of anisotropic 
magnetic excitations in the k-space.©@ These compounds are also found to be very sensitive 
to a small amount of impurities such that the residual DOS at zero energy excitations increases 
drastically in roughly proportional to the square root of impurity concentration. 

In this paper we develope a theory of k-space approach to understand the anomalous properties of 
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anisotropic semiconductor of heavy fermions, such as CeNiSn. We follow the formalism developed 
for the Fermi liquid theory of heavy fermions on the basis of the periodic Anderson model, while 
we apply it to the filling corresponding to the band insulator. An essential point of our theory 
is that the hybridization matrix element can happen to vanish along some symmetry axis (in the 
k-space) of the crystal for a particular symmetry of the crystal field which is expected to realize 
in CeNiSn. Then, the hybridization gap also vanishes along the same direction, which can explain 
wide range of the anomalies observed in CeNiSn. 

An idea that the node of the hybridization can possibly give rise to the pseudo-gap structure 
of quasiparticles has been addressed by Kagan, Kikoin and Prokof'ev,© although their theory 
seems not to have been fully developed and includes some ambiguity. For instance, the condition 
for the appearance of the node of the hybridization was neither specified, nor the explicit form of 
the quasiparticle DOS was given. Hereafter we develop the theory on a more solid ground of the 
formalism. 



We develope the formalism of our theory in §2.1 ~ §[2.6|, and discuss about DOS in §2.7 and the 



effect of impurity scattering on DOS in § |2.8| . §|j| is devoted to the discussions of physical proper- 
ties and those validity: the specific heat (§ |3.1| ), the longitudinal NMR relaxation rate (§ |3.2j ), the 
magnetic properties (§ |3.3| and § [3.4|) , the neutron scattering intensity (§|3.5|), and the anisotropic 
temperature dependence of the resistivity (§ |3.6[ ) are discussed in detail on the model calculations. 
Futhermore, the effect of pressure (§[Oj) and the lifetime of quasiparticles (§|3.8|) are briefly dis- 
cussed. 

§2. Theory 

2. 1 Hamiltonian 

We start with the periodic Anderson model keeping it in mind that (4/) 1 configuration is realized 
in Ce 3+ ion in those compounds® 

H = H c + H f + H c _f, (2.1) 

where H c , Hf, and H c -f stands for the Hamiltonian of conduction electrons, /-electrons, and the 
hybridization among them, respectively. The first term in eq. Q2.1| ), H c , is given by 

Hc = J2^ c L c ^ ( 2 - 2 ) 

where (c^ a ) creates (annihilates) a conduction electron in a plane wave state labeled by wave 
vector k and spin <r(±). The plane wave state can be expanded around site i (r^) as follows: 

1 A 00 ' 

|ka) = -/=e ikr X(T = ^e ik - r ^i^(fc|r-r,|) £ 1^* OW^r-rJx., (2-3) 
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where Xa is the spin function, ji(kr) is the spherical Bessel function, Y™ is the spherical harmonics 
with the argument of solid angle f2 r of the position vector r or of the wave vector k, and V is 
the volume of the crystal. 
The second term in eq.([2.1|), Hf, is given by 

H f = ^2 E Mft M fkM + 2^ E fiAlfiMfiM'fiM', (2-4) 
kM i,M^M' 

where f} M (fiAi) is the creation (annihilation) operator of the /-electron on the orbital M in the 
4/ shell at site i, and /j[ M (/kM) is its Fourier transform, and Em's denote the energy levels of the 
4/-electron which are split into j=7/2 and j'=5/2 multiplets under the spin-orbit interaction and 
further separated by the crystal-field interaction in general. The angular part of its eigen function 
is expressed as 

\M) = E b ff E aL^ m (a- r Jx CT , (2.5) 

/i ma 

where [i is the z-component of the total angular momentum, j, and a^ m(y are the Clebsch-Gordan 
coefficients, and b^f are coefficients specifying the crystal-field level. The last term in eq.(2.4) 



represents the Coulomb repulsion U between /-electrons in the states \M) and \M). For simplicity, 
we neglect M- and M'- dependence of U . 



The last term in eq.(2.1), H c _f, describes the hybridization between /- and conduction electrons: 

H c-f = E (V^mAJ^m + h.c), (2.6) 
kA/o- 

where V^ Ma is the mixing matrix element which can be calculated with the use of eqs.( |2.3D and 

(El) as 

V kM „ = ^V kl E bff E «L^r(^k)- (2-7) 

fjj m 

Here, V^i has the energy dependence of the mixing matrix, and is treated as a parameter of our 
model. 

2.2 Hybridization and Green Function 

Now we consider the case in which the crystal-field splitting is so large that we can neglect the 
effects of excited crystal-field states in the relevant low-temperature and low-energy phenomena. 
Then we are left with two conduction bands (±cr) and two /-levels(±M). In this case the Green 
functions of conduction- and /-electrons are given by^ 
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u — Ef — £/(k, u) 



if — 1 



G J M (M,u) = (2.8b) 



u> - , _ , , - , 

^-Ck 
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where E* is the lowest crystal-field level and S^(k, ui) is the self-energy of /-electrons due to the 
Coulomb repulsion U. The hybridization V^ 2 (k) can be regarded as independent of a after an 
appropriate linear combinations of \M) and \M) have been taken. 

Vf (k) = |Vk M(T | 2 + \V^ la \ 2 = |Vk M(T | 2 + \V kMa | 2 (2.9) 

The k-dependence of Vj(k) reflects the symmetry of the lowest crystal- field level in general.il' 0) 
For example, in the case of CeNiSn and its isostructual compounds, in which the approximate 
local symmetry of the crystal field is trigonal /^J^, the energy levels split into three doublets: 
|5/2, ±3/2) ,o|5/2,±l/2) + b |5/2, =F5/2), and -6|5/2,±l/2) + a |5/2, =f5/2), with a and 6 being 
appropriate constants satisfying a 2 + 6 2 = 1. Depending on the symmetry of the crystal- field 
ground state, there occurs various angular dependence of the hybridization V^(k). While the 
hybridization is finite at any direction of k = k/|k| in the state including 1 5/2, ±1/2) , such as 
|±p) = a |5/2, ±1/2) + b 1 5/2, =p5/2), it vanishes along the quantization axis (z-axis) in |±m) = 
1 5/2, ±3/2) 

Vgp(k) = V 2 [a 2 2(5^ - 2A; 2 + 1) + 6 2 5(1 - k 2 z ) 2 - 4 v / 10a6(A: 2 - 3^ 2 )^4] , (2.10a) 

y ± 2 m (k) = F 2 (l - fc 2 )(l + 15fc 2 ), (2.10b) 

where V 2 gives |k|-dependence of ^ 2 (k) and z-axis is taken as parallel to the o-axis, the symmetry 
axis of these crystals. 

2. 3 Effective Hamiltonian for Quasiparticles 

We are interested in the low temperature region, in which the physical properties can be described 
by the renormalized quasiparticles near the Fermi level after the many-body effect due to the on- 
site repulsion U in eq. fl2.4D has been taken into account. These quasiparticles are described by the 
effective Hamiltonian 

Heff = E £kcL c k. + £ E f flf k ± E E ( Wi/k + h.c), (2.11a) 



k<r M=± 



where 



% = ^ k [£7 / + E / (k,0)], (2.11b) 

VkMtT = VzkYkMa, (2.11c) 

where the renormalization amplitude % is defined as 



1 dE f (k,u) 
duo 



-i 



< 1 (2.12) 

LU = 



Here, the renormalized /-level Ef has a k-dependence through that of the self-energy S/(k, 0) in 
general. However, we first investigate the case where the k-dependence can be neglected, because 
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the heavy quasiparticles themselves would not be formed if Sy(k,0) had appreciable dispersion. 
The effect of its k-dependence will be discussed later in relation to the resistivity (§ |3,6| ) and the 
effect of pressure (§ |3.7| ). 
Then we can rewrite the Green functions, eqs, (|2.8| ), as 

1 



G£(k }W ) 



G f M (k,u) 



^;(k) | ^i(k) 

w-&-V/(k)/(w-£>) u-E+ u-E£' 



u-Ef-Vf 

where E^ are two hybridized-quasiparticle bands given by 

1 



4(k) At(k) 



E S-2 



t k + E f ±j(t k -E f )' + AV?(k) 



where 



V 2 (k) = z k V 2 (k). 
The residues A^_(k) and -A±(k) in eqs.( 2.13j ) are 



A c ± (k) 



4(k) 



1 + 



y f 2 (k ) 



i + 



f/(k) 



-i 



(2.13a) 
(2.13b) 

(2.14) 
(2.15) 

(2.16a) 
(2.16b) 



(^-e k ) 2 . 

which give the spectral weight of conduction- and /-electrons, respectively, in the upper /lower 
bands. 

2.4 Band Insulator of Quasiparticles 

In the case of the electron filling corresponding to the band insulator, the renormalized Fermi 
level is located in between E£ and forming the hybridization gap. As we have mentioned above, 
almost all the compounds called "Kondo insulator" contain even number of electrons in the unit 
cell and have a right to be a band insulator. A difference from the conventional semiconductor is 
that the hybridization gap is highly renormalized by strong correlation effect between /-electrons. 
Hereafter, we investigate the case where the ground crystal-field level is 1 5/2, ±3/2), which turn 
out to be consistent with anomalous properties of CeNiSn as discussed below.ll* Then, due to 
eqs.(2.10b),( 2.14j ) and ( 2.15) ) , the hybridization gap A(ke) is given by 



A(k B ) * 2^-^ ~ 2z kB ^-^ = T K (1 - k 2 Bz )(l + 15k 2 Bz 



(2.17) 



where k B denotes the wavevector at the zone boundary and Tk = 2z^ B V 2 / D. For simplicity we 
neglect |k|-dependence of V 2 and k-dependence of z^ B . In deriving eq.( 2.17 ), we have assumed that 
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the renormalized hybridization ^fz\~V is much smaller than the bare band-width of conduction 
electrons 2D. Thus the hybridization gap vanishes at points on the zone boundary where k-Qz = ±1 
and becomes a pseudogap. That is, the renormalized DOS have no clear gap threshold. This is to 
be compared to the "axial-like gap" in anisotropic superconductors, although the resultant DOS is 
totally different as discussed below. 

2.5 Mass Enhancement Factor 

The renormalization amplitude z for the particle-hole symmetric case has been calculated by 
Rice-Ueda,!) on the basis of the Gutzwiller approximation, and by Shiba0 , on the basis of vari- 
ational Monte Carlo calculations for the Gutzwiller ansatz, with the use of a model hybridization 
V 2 (k) = V 2 . A similar but more extended result has recently been obtained by numerical renor- 
malization group method®, quantum Monte Carlo calculation© and the exact diagonalization 
method© in d = oo system. We have performed the calculation similar to Rice-Ueda's with 
anisotropic hybridization U 2 (k) = V 2 (l — k 2 ), a simplified version of ( 2.10b| ). The result for the 
filling corresponding to band insulator is 

19/12^2 , 3UD s 

■ — ex P(-^) ( 2 - 18a ) 



W z \ 32V 

which is compared with that of Rice-Ueda 



D 2 ( UD \ , 

* = ^ ex H~8yV (2 - 18b) 

where the hybridization gap is constant and fully opened. 

In those model calculations, the particle- hole symmetry is assumed, so that the occupation num- 
ber of /-electron nf per site is exactly unity, i.e. nj = 1. However, this constraint is easily relaxed 
by introducing the asymmetry of conduction band on the position of the /-level. Therefore, it 
is possible to calculate the mass enhancement factor in the way similar to above not only in the 
Kondo regime, where rif — 1, but also in the valence- fluctuation regime. 

2. 6 Effect of Impurity Scattering 

It can be shown, on the basis of the Ward identity arguments, that the s-wave impurity potential 
u is renormalized by many-body vertex correction as9 

- as/(k,(j )i i 



u — > u = u ■ 



u (2.19) 

oj=0 z 



This renormalization is shown in Fig. 1 in terms of the Feynman diagram. Then, for strongly 
correlated systems where z~ l = m*/m 3> 1, the impurity scattering always becomes that of 
unitarity limit, i.e., uNp 3> 1, even if the bare potential u is moderate one, i.e., uN-p < 1. Then 
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the Green functions of the conduction- and /-electrons are given by 



1 



V 2 (k) 



G 



04 — Ef — iz k Im£ imp (u;) 
~k 



£/ - iz k ImS imp (a;) 



V> 2 (k) 



(2.20a) 



(2.20b) 



where the self-energy Si mp (o;) due to impurity scattering is given in the t-matrix approximation by 

u 



J imp 



{04) 



imp 



1 (k 



(2.21) 



04 



where nj mp denotes the impurity concentration. In deriving eqs.( |2.20] ), the self-energy of conduc- 
tion electrons due to impurity scattering has been neglected because the renormalization, such as 
eq.( 2.19| ), does not occur. Equations ( 2.20| ) and ( 2.2 1] ) need to be solved self-consistently as in 
the case of heavy fermion superconductors, where the impurity scattering in the unitarity limit is 
known to give rise to appreciable residual DOS in the V-shaped gap even for a very small impurity 
concentration .H' 

2. 7 Density of States of Quasiparticles 
The quasiparticle DOS are calculated as follows: 



k. 
N F 



dk, 



D 



dE\l + 



D 



vf{k 

E 2 



6(o4-E)e(\E\-A(k z )), 



where 



A(k z ) = T K (l-k 2 z )(l + l5k 2 z ). 



(2.22a) 
(2.22b) 

(2.23) 



In deriving eq. (2.22b) from eq. (2.22a), we have assumed for simplicity that the conduction band 
has a linear dispersion with constant DOS, Np, and extending from —D to D, and Et = 0. A 
result of numerical calculation of N(uj), eq.(2.22b), is shown in Fig. 2. Here, the relation between 
the hybridization V and renormalization factor z is determined as zV 2 /D = 0.01D. The shape of 
DOS shown in Fig. 2 has two characteristics: (1) N(lo = 0) is finite, and (2) it exhibits four-peak 
structure, i.e., there exist two energy scales (Ax = 7k = 0.02D,A2 — 0.08D). The former (1) 
results from the fact that the hybridization gap vanishes at points {k z = ±1). The shape of DOS 
around the Fermi level 04 = can be calculated analytically as 



N(uj) ~ N F 



D 2 



64zV 2 



1 + 2 1 + 



19D 2 



12 • 32zV 2 / D 



04 

— + 



(2.24) 
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This is in marked contrast with the case of heavy-fermion superconductors, where the point node 
leads to DOS proportional to lo 2 . The reason for the characteristic (1) to hold is that there exists 
a singularity oc E~ 2 , in the first factor of the integrand of (2.22b), which arises from the Jacobian 
|d£/d£| = l + V 2 (k z )/E 2 . 

The latter characteristic (2) is related to the existence of two extremum values of A(fc 2 ), eq.( |2.23j ): 
Ai corresponds to the minimum of eq.( p.23| ) at k z = 0, at which the hybridization gap is given by 
T%, and A2 corresponds to the maximum of eq.( f2.23| ) at k z = y/7/15, in which the hybridization 
gap becomes maximum. The ratio of A2 and Aj is given as A2/A1 = 64/15. 

2. 8 Effect of Impurity Scattering on Density of States 

Next let us consider the effects of impurity scattering on quasiparticle DOS. Prom eqs.(2.20b) 
and (|2.21|), DOS is calculated self-consistently as follows: 




-z k ImS imp (tj) 



z k V 2 (k) Y 
Ef J —— + [z k Im£ imp (w)J 



(2.25a) 



ImE imp (u;) = -n imp tl^M (2 .25b) 



1 + U 2 [tTZ^N^J 



Results of numerical solution of eqs.(2.25) are shown in Fig. 3. One can see that the residual DOS at 
the Fermi level, N(tu = 0), is very sensitive to the impurity concentration and drastically increases 
with the impurity concentration ni mp . In Fig. 4 one finds that N(co = 0) is roughly proportional to 
v^frnp. Precisely speaking, the residual DOS exists, even if n; mp = 0. So, in the limit, n\ mp — > 0, 
its nj mp -dependence is given by + con; mp , where ^0 is the residual DOS without impurities, 

and Co is a proper constant of order unity. This is to be compared with the impurity effects in 
heavy fermion superconductors .011^ The square-root dependence of the residual DOS has also 
been derived on the basis of a different picture, where it is understood as an impurity band similar 
to the doped semiconductor; namely the doping accompanied by variation of carrier number is 
necessary to obtain finite DOS in the true gap. Our theory has been developed to discuss the 
case where carrier number does not change, while it is easily extended to the case where carriers 
are doped. 

§3. Physical Properties 

In this section, we study the qualitative aspects of several physical quantities and compare them 
with experiments. 
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3.1 Specific Heat 

The specific heat is calculated on the basis of the quasiparticle picture, and the electronic specific 
heat coefficient, 7 = C(T)/T, is given as follows: 



00 



7 = 2/ dxN(E)x 2 sech 2 (x), (3.1) 







where x = E/2T. The temperature dependence of 7 is calculated with the use of eqs. (fT22|) and 



( p. 25 ) and is shown in Fig. 5. A peak structure is found at T ~ Ai/2. It is noted that 7 is finite 



at T = 0, which results from the existence of the residual DOS at n; mp = 0. If Ai is fixed as 
Ai/2 = 7K, these results are in good agreement with the experimental dataU'Elil' at T < A\, 
where the maximum of theoretical curve for 7 is adjusted so as to agree with the experiment of 



ref. 23. We have also verified that the same quality of agreement with the data of ref. (24 is obtained 



while its absolute value of C m /T is about 10% larger than that of ref. 23. However at T > Ai 
the agreement become poor. It is partly improved by considering temperature dependence of the 
quasiparticle DOS. 

Futhermore, assuming that the effect of the magnetic field is only inducing the Zeeman splitting, 
the magnetic-field dependence of 7 is obtained as 



/ dxN(E) x 2 + sech 2 (x 2 + ) + x 2 _sech 2 (x 2 _) , (3.2) 
Jo 1 J 



where x± = (E± h)/2T, h = gj[iB\Jz\H , gj being a g-factor. The results for various temperatures 
are shown in Fig. 6. The coefficient 7 at low temperature exhibits two-peak structure, which reflects 
the peak structure of DOS. This prediction has not yet been observed, partly because the strength 
of the magnetic field is not enough. 

3. 2 Longitudinal Relaxation Rate of NMR 
The longitudinal NMR relaxation rate, 1/Ti, is obtained as follows 

J-oclimV 1111 ^^ , (3.3a) 
TIT u-*o^ u K J 

roo TP 1 

oc j o d EN(E) 2 sech\—)-. (3.3b) 

Here we have assumed that the quasiparticle DOS directly affects 1/T\ at Sn site via the c-f 
exchange as done in the analysis of experimental data. The results of numerical calculations of 
(3.3b), together with experimental data,© are shown in Fig. 7 for the same parameters as in Fig. 5. 
For the temperature region, O.lAi < T < Ai ~ 0.02D, 1/Xi shows the T 3 -like behavior reflecting 
the formation of the pseudogap at u ~ Ai/2 of DOS as shown in Fig. 2, and for T < 0.1 Ax, it shows 
T-linear behavior reflecting the residual DOS at the Fermi level. These behaviors well reproduce 
the T-dependence of 1/Ti observed in the experiments,© if Ai is fixed as Ai = 14K. In addition, 
as increasing impurity concentration n; mp , the residual DOS rises up drastically and the T-linear 



10 



behavior masks the T 3 -like behavior. These are also in agreement with the experiments,!^ where 
the residual DOS shows the ni mp -dependence quite similar to the theoretical curve shown in Fig. 4. 

The above results have been derived on the basis of the quasiparticle picture, so that, strictly 
speaking, its validity is assured only in the low temperature region T < Tk. However, it may be 
extended to much higher temperature region as far as the qualitative aspects are concerned. Indeed, 
1/Ti in Fig. 7 exhibits the localized character of /-electrons for T > Tk = 0.02D, 1/Ti oc T°, and 
approaches asymptotically to the Korringa-like behavior at much higher temperatures, \jT\ oc T, 
as can be inferred from the DOS of quasiparticles shown in Fig. 2. In the latter region, the NMR 
relaxation is expected to occur mainly through the coupling with the conduction electrons as in 
LaNiSn.© Recently such a behavior has been recognized by an analysis of the data of 1 /T\ up to 
the room temperature.© 

3. 3 Magnetization 
The quasiparticle contribution to the magnetization is calculated as 

N(E) 



M 



r dE [f(E -h)- f(E + h)] -L±. (3.4) 

J — oo ^ 



The results are shown in Fig. 8 where one can see that the magnetization, in the unit qj^b\Jz\, is 
proportional to h, in the unit D, at low magnetic field and the slope is given by the residual DOS. 
However, as increasing h the magnetization drastically increases at h > 0.01D owing to the two 
large humps of DOS at uj = ±Ai/2. These tendencies are found in the experimental data©. The 
slope of the magnetization at low field is enhanced by small amount of impurity and the whole 
structures of M-h curve shade off. 

3.4 Uniform Spin Susceptibility 

The uniform spin susceptibility along the easily axis (a-axis), in the unit {gj^,-Q\J z \) 2 / D, is given 
by the derivative of the magnetization as 

dM 



Re X (0,0) 



dh 



poo 

= / dE 


[ df(E)] 


h^O J-00 


dE 



N(E). (3.5) 



Its temperature dependence shown in Fig. 9 exhibits the peak structure like 7 as discussed in §3.1 
However, the temperature at which Rex(0, 0) has the maximum (T ~ Ai) is different from that for 
7. Increasing impurity concentration, the sharp dip at low temperature is filled up rapidly. Since 
the susceptibility (^^) is given only by the contribution of the quasiparticles, the Van Vleck term 
is not included. If the latter is simply a constant, the observed Knight shift@) for H\\a represents 
the behavior of the uniform spin susceptibility x a > which is in agreement with the curves of Fig. 9. 
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3.5 Neutron Scattering 

The imaginary part of the dynamical susceptibility, the spectral weight of spin fluctuations, is 
calculated at T = without a vertex correction as 

Im X (Q, oj) ~ tt^ [f(E k ) - f(E+ +Q )} S(oj - E+ +Q + E£) (3.6) 
k 

where f(E) is the Fermi distribution function and f(E^) = 1 and f(E^) = at T — > 0. The 
spectral weight (|3.6|) has been calculated numerically with the use of the quasiparticle dispersion, 



eq.(2.14), for specified Q's. The spectral weight at Q = (1/2,0,0) is shown in Fig. 10 (a), which 
shows that there exists a broad hump at around oj = 0.08-D. (It is noted that x-, y-, and z-axis 
here corresponds to the b-, c-, and a-axis, respectively of CeNiSn; so that Q = (1/2,0,0) implies 
Q = [0, 1/2,0] in the notation of experiments of CeNiSn for instance.) This structure corresponds 
to the transition from one peak of DOS at oj = — 0.04D to another at oj = 0.04D in DOS. For 
example, the former peak corresponds to the quasiparticle at ki = (—1/4,1/2,1/2), while the 



latter at k 2 = (1/4,1/2,1/2), because k Xz = k 2z ~ VT/TE. Thus, Q = k 2 - ki = (1/2,0,0). 
The spectrum at Q = (0,0,1/2) shown in Fig. 10 (b) has a peak at around oj = 0.06D, which 
corresponds to the energy from an edge of the gap at oj = — 0.02D to one peak at oj = 0.04D in 
DOS. For example, the former peak corresponds to the quasiparticle at k3 = (1/2,0,0), while the 



latter at k 4 = (1/2,0, 1/2), because k 3z = and k 4z ~ y/7/15. Thus, Q = k 4 - k 3 = (0,0, 1/2). 
It is seen that a difference between the two spectra arises from that of the possibility of the zero 
energy transition. Since the hybridization gap vanishes along the z(a)-axis, there exist zero energy 
excitations for the transition between the points on the z(a)-axis. This is because we make choice 
of non-dispersive /-level. However, generally speaking, Imx(Q,0) have to vanish in the normal 
Fermi liquid theory by the symmetry reason. This anomaly is removed by taking into account the 
small dispersion of /-electrons due to a possible weak k-dependence of the /-electron self energy, 
S/(k,0). 

The spectral intensity at oj = 0.08-D ~ A 2 and Q = (Q x ,0,0) = [0, Qb, 0] is computed as a 
function of Q x (Qb) and shown in Fig. 11. One can see the peak at Q x (Qb) = 1/2- The reason is 
that the peak shifts to higher energy as deviating from Q x (Qb) = 1/2. These features are consistent 
with Qfe-dependence of the intensity at 4.25meV in the inelastic neutron scattering.© 

The details of these spectra are modified according to choices of the dispersions of conduction 
electrons, though characteristic structures do not change. However, when the conduction band do 
not cross the /-level in some direction in the k-space, the excitation energy is rather higher and 
the spectrum can not be observed at low energy region in general. This case may be realized in 
the low-energy spectrum at Q || c-axis in the experiment.® 
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3.6 Anisotropy of Resistivity 

It is the resistivity that is one of the measures to classify the heavy- fermion materials into "Kondo 
insulator" or not. The resistivity in heavy fermions exhibits the Kondo effect at high temperature 
region and metallic or activation-type behavior at low temperature region. We are interested in 
the behavior at temperatures lower than the coherent temperature T co h, in which the current is 
carried by the quasiparticles. In this case the conductivity can be evaluated by 

<V(T) oc £ J kM J k ,r k ("^^) , (3-7) 

where J k/J is the velocity of the quasiparticle and r k is its lifetime.© Assuming that the microscopic 
expression of the current is given only by the conduction electrons (neglecting the dispersion of 
/-electrons), the anisotropy of the conductivity is given as follows: 



k V dE k 



(3.8a) 



x / ,lE k / dn k ^ , ^ f- dffi" j ■ (3.8b) 

(-z k ImS / (k,ii; k ))y / 2 (k) + (-ImS c (k, J E k ) J B2) V dE k 



( T ± OC^^(k)V(k) 2 7 k (- 

k V 



(3.8c) 



J J f-a.ImS f fk. E,)) V?(\c) + -\m^Jk.E,,)E?) \ dE k J 



(-ZylmXfiKE^VfQi) + (-ImScCk,^)^ 

where cry and a± are the conductivity along the a-axis and in the 6c-plane, respectively, and 
1 1 



Tk El + y?(k) 



(-^ImS^k,^)) V> 2 (k) + (-Im£ c (k,£ k )£2)j ; ( 3 . 9a ) 

Vz{k) = W z and ^' (k) = w x - (3 - 9b) 

Here we have introduced E c (k, £" k ), the self-energy of the conduction electrons, because a» in 
pure system diverges logarithmically otherwise, reflecting the fact that the conduction electrons 
are decoupled from /-electrons in the z-direction (a-axis) where the hybridization vanishes. We 
regard the renormalization factor of conduction electrons from E c (k, E k ) as 1. In this model the 
conductivity in the 6c-plane is isotropic unless the anisotropy of the conduction band is taken into 
account. In deriving (3.8d) from (3.8c), we have taken into account only the low energy excitations 
around k z = ±1 so that the obtained result should be regarded as that for asymptotic behavior in 
the limit T — > 0. 

In order to discuss the temperature dependence of the conductivity at low temperature region, 
we must calculate the energy dependence of the imaginary part of the self-energies. For simpicity 
we calculate these along a standard treatment of the Fermi liquid theory neglecting the momentum 
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dependence of the full vertex: 

Im£ M (k, E+) oc J dpdq^(k - q)^(p)^(p + q)6(E+ + E~ - £+_ q - £+ +q ), (3.10) 

where \i = c or /. These integrations are computed by the Monte Carlo calculation. The numerical 
results are shown in Fig. 12. From this one can see that the energy dependence of ImS c (k, EjJ 
and ImSj(k, EjJ near the Fermi level at zero temperature can be approximated by E^ and E^, 
respectively. Futhermore, we assume that the temperature dependence of ImS«(k, EjJ is given with 
replacing E^ by max (E?_ and (7rT) 2 ), as can be seen from the structure of the Green functions. 
By using these results we can estimate the temperature dependence of the resistivity. Substituting 
these energy dependence into eqs. fl3.8|) , we obtain up to the logarithmic accuracy 



an oc T 1 , and a± oc T°, (3.11a) 

i.e, 

pn oc T , and p± oc T°. (3.11b) 

Next we discuss the effect of impurity scattering on For simplicity we first assume that the 
self-energy Im£ c (k, E^) and ImSj(k, E\f) are independent of E^ and proportional to the impurity 
concentration n mip . In this case, by using eq.(3.8b), we obtain 

crj| mp oc T 2 /n imp , i.e. p™ p oc n imp T" 2 . (3.12) 

However, if the small dispersion of /-electrons due to a possible weak k-dependence of the /-electron 
self energy, £f(k, 0), the current can be carried also by /-electrons, so that 



Cm oc 



dEk 2^ (Jjm\ , (3 . 13 a) 

J ^(-z k ImS / (k,S k ))V> 2 (k) + (-ImS c (k,ii; k K) V 0E k J 

oc . (3.13b) 

Thus the singularity of the resistivity (|3.12j ) at T = is suppressed. Nevertheless, the residual 
resistivity in the limit T — > increases drastically as increasing the impurity concentration Ti[ m p. 
If we take all the contributions of the quasiparticles into consideration, including the logarithmic 
corrections, the resistivity is given by 

"ii = ; i—fi r=r. P.") 

K + PW) « + (,„ g ( C2rK K /r) + (r/T K )C( Cs rK/r)j 

where ci(~ 0.1), C2(~ 10) and C3(~ 5) are fitting parameters, which are connected with the small 
dispersion of /-electrons, the interaction between conduction electrons and impurity scattering of 
conduction electrons with the Born approximation, respectively. ni mp /3o/ci is the resistivity at 
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T — > 0. The resistivity for proper parameters are shown in Fig. 13. These results are in qualitative 
agreement with the experimental dataS© In particular, the temperature dependence of p\\{p a ) 
observed in the best sample to date is well reproduced as seen in Fig. 13. It is also found that the 
resistivity is sensitive to the concentration of impurities at low temperature in consistent with the 
experiments. 

It should be remarked here that the weak k-dependence of E/(k, 0) inevitably gives rise to a 
small semimetallic Fermi surface around X-point (intersection of a-axis and the zone boundary) 
in general, so long as the hybridization vanishes along the a-axis as eq.(2.10b). However, such 
a small Fermi surface is expected to give only little effect on the qualitative behavior of DOS 
of quasiparticles discussed in previous sections, while it sensitively affects the low temperature 
behavior of the resistivity especially in the case where the impurity scattering greatly enhances the 
resistivity as in eq. ( |3.12 ) when there exists no k-dependence of S/(k, 0). In deriving ( 3.14| ), we have 
taken into account the dispersion of /-electron through the k-dependence of £/(k, 0), nevertheless 
we have used the same DOS as eq. ( |jg ). In this sense, the calculation is not self-consistent and 
the resistivity ( 3,14| ) should be regarded as a provisional one. However, the expression ( |3.14j ) gives 
a good description for T ~ Ai/2, or a good starting point at least. 

We should have a few words about p±. The resistivity pb and p c of the same sample as shown 
in Fig. 13 exhibits a dip at T ~ 3K and saturate at T —* 0. This may be understood as follows. 
Since quasiparticles around k z = with excitation energy (~ Ai/2) have large dispersion along 
x- or y-direstion, those are expected to contribute considerably to the conduction perpendicular 
to z(a)-axis at T ~ Ai/2, leading to the suppression of pj_ (pt, and p c ). If we use Ai/2 ~ 7K 
estimated above by means of 1/Ti's result, such suppression or dip is expected to occur at T < 7K 
in consistent with the above observation. 

Now we briefly discuss about the longitudinal magnetoresistance (MR). The change from the 
positive MR to the negative MR is observed with increasing mag netic fieldJUi) The negative 
MR at higher field region may be caused by the suppression of spin fluctuations by magnetic field 
such as in the impurity Kondo effect. However, the case of the positive MR at low field is more 
complicated. We have calculated magnetic-field dependence of the lifetime t\, of quasiparticles and 
verified that is a decreasing function of the magnetic field in the low-field region. Although 
DOS at low energies increases with the magnetic field due to the Zeeman splitting of DOS, shown 
in Fig. 2, the weight of conduction electrons of those states decreases in general leading to the 
enhancement of the resistivity. Therefore, MR is determined on such a delicate balance between 
the effects on the lifetime of the quasiparticles and the details of DOS at the Fermi level. The 
experimental data can be understood as follows: MR at lower fields is positive by shortening of the 
lifetime and at higher fields becomes negative due to the drastic increment of quasiparticles which 
can carry the current. 
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3.7 Pressure Dependence 

Next we discuss the pressure dependence on the property of quasiparticles. As the lattice con- 
stant becomes short under the pressure, both the band-width D of conduction electrons and the 
hybridization V are enlarged. However, the fundamental energy scale V 2 /D is expected to be an 
increasing function of the pressure because V is much more sensitive than D for heavy fermions 
where V arises through rather small overlap between /- and conduction electrons. Much more 



pronounced effect of pressure on the hybridization gap (2.17) arises through enlargement of the 
renormalization amplitude zy,. This is because % is an exponentially small qauntity as (2.18) for 
heavy fermions so that its relative change under pressure is far larger than that for V and D them- 
selves. Therefore, the energy scale of the gap is expected to increase by applying the pressure, and 



so is the resistivity in eq.(3.14), which is scaled by 2k • 

This tendency is consistent with the behavior of CeNiSn and CeRhSb where the peak of the 
resistivity shifts to the high temperature by the pressure. 00 And also the suppression of those 
resistivity with the pressure in the limit T —* can be understood as follows. As discussed in § |3.6| , 
there exists a very tiny semimetallic Fermi surface at around k = (0, 0, ±1) in general, so long as 
£/(k, 0) has the dispersion along the a-axis no matter how small it is. After the Fermi surface 
grows further under the pressure due to the growth of the dispersion of S/(k, 0), an apparent 
semimetallic behavior is expected to prevail leading to the suppression of the resistivity. That is to 
say, a parameter ci, which corresponds to <9Xy(k, 0)/<9k, increases with the pressure in the formula 



( p. 14 ), so that the resistivity in the limit T — ► 0, ni mp po/ci, is suppressed. 



3.8 Lifetime of Quasiparticles 

We have neglected an effect of quasiparticle damping due to inelastic scattering so far. Here we 
briefly discuss its effect on DOS of quasiparticles and temperature dependence of physical quantities. 

According to eq. d3.9a| ), the lifetime of quasiparticles around k z = ±1 is nearly proportional to 
(E^ + VfOt)) /Vf(k)E^ at low energy region E^ < A\. This lifetime is longer than in the case of the 
normal Fermi liquid theory for E^ < Ai. This is because the scattering between quasiparticles is 
suppressed at low energy, owing to the restriction of phase space satisfying the energy-momentum 
conservation law. Namely, the quasiparticles with low energy are located around (0,0, ±1), so that 
such phase space is restricted within narrow region around (0,0, ±1). These quasiparticles make 
the flat part near the Fermi level at DOS of Fig. 2. 

The quasiparticle, corresponding to the peak structure at u = Ai/2 in DOS of Fig. 2, are located 
along k z = 0. These quasiparticles also suffers little inelastic scatterings, again because of the 
restriction due to the energy-momentum conservation law. Thus these lifetime is very long, leading 
to ImT, f(k z ~ 0, E\g) ~ 0. So it is expected that the peak structure at u = Ai/2 in DOS remains 
even if the effect of inelastic scattering is taken into account. 

On the contrary, we have no reason to keep the sharp peak structure at to = A2/2 in DOS, 
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because the restiction due to the energy-momentum conservation does not suppress the inelstic 
scattering of quasiparticles forming this peak. So the structure at u = A2/2 probably becomes a 
broad hump. 

However, we believe that the two peak structure in DOS of Fig. 2 remains even though the 
effect of inelastic scattering is taken into account. Indeed, the DOS calculated by the second order 
perturbation theory exhibits two peak structure similar to those of Fig. 2, although the peak at 
lo = A2/2 is somewhat broadened. This result will be discussed elsewhere. 

Futhermore, T > is the temperature region where the damping effect of the quasiparticles 
affects temperature dependence of physical quantities. For T > Ai, these peak structures of 
DOS fade out, while the physical quantities are averaged by temperature dependence of the Fermi 
distribution. Therefore, the neglect of the damping effect may give rise to no serious errors as far 
as the qualitative temperature dependence is concerned. 

§4. Summary and Discussions 

On the basis of the idea of "adiabatic continuity", a theory of the anisotropic semiconductor 
of heavy fermions has been developed to explain the anomalous properties of CeNiSn and its 
isostructural compounds. A difference from the conventional semiconductors is that the band gap 
is formed by the highly renormalized quasiparticles near the Fermi level. So the gap has meaning 
only at low temperature region T < Ai (corresponding to the hybridization gap), while the coherent 
peak of quasiparticles fades out exhibiting the behaviors of the Kondo lattice metals.© 

Wide range of anomalies of CeNiSn can be understood by a model of the anisotropic hybridiza- 
tion gap which vanishes along the a-axis. The anisotropy of the gap reflects a k-dependence of 
hybridization matrix elements between the conduction electrons and the /-electron with particular 
symmetry of the crystal field state. The desirable k-dependence occurs if the lowest crystal field 
state consists mainly of 1 5/2, ±3/2) due to its approximately trigonal symmetry and the conduction 



electrons near the Fermi level are described by the plane waves, as discussed in §[2T] and § |2.2 . 

Since there is no band calculation of LaNiSn available to date, it is difficult to assess whether the 
latter condition is fulfilled in CeNiSn. However, it may be not unrealistic to assume that the state 
of conduction electrons hybridizing with the /-electron localized at Ce site can be approximated by 
the plane waves with the wave vector |k+G| < several x (2n/a), G being some reciprocal lattice 
vector and a being the lattice constant. This is because the only way for the /-electron to mix with 
electrons on different sites is through the mixing with the plane wave states outside the muffin-tin 
spheres so long as the conventional LAPW calculation is performed. 

Band calculations of CeNiSn shows that the bands around the Fermi level have mainly Ce 4/ 
character with mixture of Ni 3d component. So, in the tight-binding picture, the hybridization 
is expected to arise through the overlap of Ce 4/ wavefunction and tails of Ni 3d wavefunction. It 
is seen by a simple calculation of the tight-binding model that the hybridization between /-electron 
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in the state 1 5/2, ±3/2) and d-electrons on the surrounding ions with trigonal symmetry vanishes 
on the k z axis, i.e., V(0, 0, k z ) = 0. 

In order to obtain more solid picutre of the k-dependence of the hybridization, we need more infor- 
mation of the band structure of LaNiSn. It is also interesting to discuss a difference between CeNiSn 
and the so-called "Kondo insulator", such as Ce 3 Bi4PtJ*l> and YbBi 2 §), with non-vanishing gap 
in any directions of the Brillouin zone. From the present point of view, its difference is attributed 
to that of the k-dependence reflecting the symmetry of the lowest crystal-field level. We leave such 
discussions for future studies. 
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Fig. 1 The diagram for the many-body vertex correction of impurity scattering. The broken line represents 
the impurity potential u of s-wave, the external solid line the Green function of the /-electrons, and internal 
solid line stands both /-electrons and conduction electrons. T is the full vertex due to the Coulomb 
repulsion between /-electrons, u is a renormalized potential of impurity scattering. 

Fig. 2 The density of states, N(u>), of the quasiparticles as a function of ui in the unit D, half the band- 
width of conduction band. The Fermi level is located at u = 0. 

Fig. 3 The effect of the impurity scattering on the density of states. 

Fig. 4 The density of states at the Fermi level as a function of the impurity concentration. 

Fig. 5 The temprature dependence of the specific heat coefficient 7. The unit of the temperature is D, half 
the band-width of conduction electrons. Circles are experimental data (C m /T) of ref. Ai/2 = 0.01-D 
corresponds to 7K. 

Fig. 6 The specific heat coefficient, 7, as a function of the magnetic field, h. The unit of h is D, half the 
band-width of conduction electrons. 

Fig. 7 The longitudinal NMR relaxation rate, 1/Ti (in arbitrary unit), as a function of the temprerature 
T in the unit D, half the band-width of conduction electrons. Triangles and crosses are experimental data 
of ref. ||. Ax = 0.02£> corresponds to 14K. 

Fig. 8 The magnetization, M (in the unit gjHB\Jz\), as a function of the magnetic field h (in the unit D) 
at T = 0. 

Fig. 9 T-dependence of the uniform susceptibility, Rex(0, 0) in the unit {g j p-e,\J z \) 2 / D . The unit of T is 
D, half the band-width of conduction electrons. 

Fig. 10 The spectral weight of spin fluctuations, Imx(Q, uj) (in arbitrary unit), as a function of lo at T = 
for (a) Q = (1/2,0,0) = [0,1/2,0] and (b) Q = (0,0,1/2) = [1/2,0,0]. Imx(Q,0) is finite in (b) due to a 
choice of non-dispersive /-level. 

Fig. 11 The intensity of Imx(Q,w = A2) in arbitrary unit as a function of Q x (Qb) of Q = (Q x ,0,0) — 
[0,Q b ,0] at T = 0. 

Fig. 12 The energy dependence of the self-energies, (a) Im£ c (k, Ej,) and (b) ImS/(k, E^) at zero tem- 
perature. The straight line in (a) and (b) shows E$ and E^ dependence, repectively. The unit of E-^ is 
D, half the band- width of conduction electrons. 

Fig. 13 The resistivity p» as a function of the t emperature T, in the unit D, half the band-width of 
conduction electrons. The parameters in eq.( 3.14 ) are chosen as po — 400, c% — 0.1,^a = 10 and C3 = 5. 
Closed circles show the temperature dependence of p a of the best sample of CeNiSn.E!rEf 
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